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We propose a new scenario for the early universe where there is a smooth transition between 
an early de Sitter-like phase and a radiation dominated era. In this model, the matter content is 
modelled by a new type of generalised Chaplygin gas [3| for the early universe, with an underlying 
scalar field description. We study the gravitational waves generated by the quantum fluctuations. 
In particular, we calculate the gravitational wave power spectrum, as it would be measured today, 
following the method of the Bogoliubov coefficients. We show that the high frequencies region of 
the spectrum depends strongly on one of the parameters of the model. On the other hand, we use 
the number of e-folds, along with the power spectra and spectral index of the scalar perturbations, 
to constrain the model observationally. 

PACS numbers: 98.80.-k, 98.80.Bp, 95.36.+X, 04.30.-w 
Keywords: Chaplygin gas, inflation, gravitational waves 



I. INTRODUCTION 

The inflationary paradigm is the most consistent 
scenario to explain the origin of the large scale structure 
(LSS) of the universe, as well as the anisotropics in 
the cosmic microwave background radiation (CMBR) 
0, Originally, an early inflationary era in the 

universe was invoked as a mechanism to solve several 
shortcomings of the big bang theory j4|. Afterwards, 
it was realised that an early inflationary era of the 
universe generates density perturbations as seeds for the 
structure we see nowadays. In addition, a background 
of stochastic primordial gravitational waves (GW) is 
also predicted and it is originated from the vacuum 
fluctuations. The latter (GW) if ever detected in the 
future would provide an amusing source of information 
about the early universe. On the other hand, and 
most importantly, the predictions of the inflationary 
theory have been corroborated by several cosmological 
observations like the recent WMAP5 data (sj- 

In this letter we propose a phenomenological model 
for the inflationary era and the subsequent radiation 
dominated epoch of the universe. In particular, we 
suggest a way to extend the framework of the generalised 
Chaplygin gas (GCG) to the first stages of the universe^. 
The GCG has attracted a lot of attention over the last 
years [l], 0413 • It was initially introduced in cosmology 
as a mean to unify the dark sectors of the universe, i.e. 
dark energy and dark matter, and has been contrasted 
with several cosmological observations ^ . On the other 
hand, it has been as well shown that some GCG models 
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^ An alternative inflationary model inspired on the Chaplygin 
gas was proposed in Ref. [3I where only the inflationary era is 
accounted for. 



can be an excellent framework to analyse dark energy 
related singularities 0, [13] • 

In this work, we investigate the possible imprints in 
the power spectrum of the stochastic background of GW, 
for the model we propose for the transition from the 
inflationary era to the radiation dominated epoch^ . This 
transition is far from being well known and it can be of 
some interest to explore the signatures associated with 
it, shedding light on the inflationary model behind the 
accelerated expansion in the early universe. 

We calculate the GW production using the Bogoliubov 
coefficients which obey a set of differential equations ^12| - 
ITsj . This method has advantages over the frequently used 
sudden transition approximation, because, associated 
with it, there is always an overproduction of gravitons 
of large frequencies, which is avoided in a natural way 
by the use of the continuous Bogoliubov coefficients [la] . 
This is also a very practical method to calculate the full 
spectrum, from the very low frequencies corresponding 
to the present cosmological horizon 10~^''llz, till those 
large, kHz up to GHz, frequencies associated with the 
transition between the inflationary and the radiation- 
dominated universe. 

The present letter is organised as follows. In section [H] 
we present our model based on a new type of generalised 
Chaplygin gas for the early universe. This fluid can 
be as well represented by a minimally coupled scalar 
fleld as shown in section IIIII In section IIVI we 
constrain our model observationally. In section |Vl we 
summarise the methodology used to obtain the spectrum 
of the stochastic GW which is based on the Bogoliubov 
coefficients. We present in section IVII the spectrum of 
the GW for the model introduced in sections HI] and IIIII 



^ In Ref. [Till , the spectrum of GW for a universe whose dark 
energy component corresponds to a Chaplygin gas was analysed. 
This scenario is different from ours. 
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Finally, in section IVlIl we present our conclusions. 

II. A MODIFIED CHAPLYGIN GAS FOR THE 
EARLY UNIVERSE 

We consider a perfect fluid whose energy density, p, 
scales with the scale factor, a, as 

^=(^+;^)"°' (2.1) 

where A, B and a are constants. This energy density 
interpolates between a positive cosmological constant at 
small scale factors and a radiation fluid at large scale 
factors as long 1 + a is negative and A, B are positive. 

By assuming the conservation of the energy density of 
such a fluid (cf. Eq. (|2.ip ). which is a direct consequence 
of the Bianchi identity, the equation of state satisfied by 
its pressure, P, and energy density, p, reads 

(2.2) 

Before moving on a few words on the physics behind such 
a fluid are in order: like a generalised Chaplygin gas, 
whose equation of state reads P = —Alp"", a fluid with 
the equation of state (|2.2p could be formally related to 
an effective description of a complex scalar field whose 
action can be written as a generalised Born-Infled action 
corresponding to a perturbed d-brane in a (d + 1, 1) 
space-time. For the propose of this work it is sufficient 
to describe such a matter content through the equation 
of state (|2.2p or through a scalar field as it is shown on 
the next section. 

The equation of state (|2.2p corresponds to a mixture 
of a radiation fluid and a generalised Chaplygin gas 
[U, 0] although, it is only the total energy density p 
which is conserved. This type of equation of state has 
been previously introduced in the context of dark energy 
models [TMl (see also [IB-il). Nevertheless, we would 
like to stress that it has never been analysed in the 
context of the early evolution of the universe. More 
precisely, the equation of state (12. 2p has been analysed in 
Ref. [3l for < a < 1; i.e. a fluid interpolating between 
a radiation at early time and a cosmological constant 
at late time. Here we will be analysing a completely 
different range of the space parameter where a < — 1. 

For a < —1 the energy density (12. 1[) can be 
approximated as 

p ~ Att^ a < a-i,, (2.3) 
p ~ a > a^, (2.4) 

where 

/ 5\ 4TTTS7 

». - (,) . (2^5) 



Therefore, a homogeneous and isotropic universe with 
this matter content (cf. Eqs. (|2.ip and (12.21) ) undergoes 
a primordial inflationary era; initially described by a 
de Sitter (dS) expansion (see Eq. (|2.3p ). Then, the 
universe keeps inflating until the scale factor reaches 
a* where inflation ceases; i.e. p + 3P = 0. Finally, 
it enters a radiation dominated epoch (see Eq. (12. 4p ). 
The transition from the inflationary era to the radiation 
dominated phase takes place smoothly. 
The Friedmann equation 

= yp, (2.6) 

where = SttG and p is deflned in Eq. (|2.ip . can be 
integrated analytically^ [2^, [13] 

_|A^<i+=) (ry - ^TF (^-r, 1; 1 - -; - j 
-y-h(x + yYF{-r^■^~-^^^L^. (2.7) 

In the previous expression r\ is the conformal time, 77* is 
a constant, F(a,b;c;y) is a hypergeometric function^ and 

. x-4(l+a) 

\a* ) 2(1 + a) 

At the conformal time ry^, the universe exits the 
inflationary era as a = a^,. 

The primordial universe starts its dS expansion at 77 — > 
—00 where a a^,. On the other hand, the universe gets 
radiation-dominated a,t ri — >■ 00 where <C o. At this 
respect, we note that [2J| 

F(-r,l;l-^;l) =-1. (2.9) 

By choosing 

^A^r;. ^ -Y F (^-r, 1; 1 - ^; ^) , (2-10) 

the relation between the conformal time, 77, and the scale 
factor given in Eq. (|2.7p can be simplified to 

Ka^A^oWTy = -VSy^^{l + yy x (2.11) 

xF('-r,l;l-^;^V 



^ This Friedmann equation can be integrated in two steps by (i) 
performing it for < 1 + o [23II and (ii) performing an analytical 
continuation of the hypergeometric function to values of a such 
that 1 + a < [23 |. 

* A hypergeometric series F(fe, c; d\ x), also called a hypergeometric 
function, converges at any value x such that < 1, whenever 
6 -I- c — d < 0. However, ifO<fe-|-c — d< 1 the series does not 
converge at x = 1. In addition, if 1 < b+c—d, the hypergeometric 
function blows up at |a;| = 1 |23|| . 
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We notice that this redefinition of 77 does not modify the 
fact that in this model the universe has an infinite past (in 
terms of the conformal time), where it is asymptoticahy 
dS, and an infinite future (in terms again of the conformal 
time), where it is radiation-dominated. 

The whole evolution of the universe can be described 
by the modified Chaplygin gas (|2.2p at early-time and by 
the ACDM model at late-time; i.e. 

{A + ^ittko ) ^ early-time 
p=^ ^ (2.12) 

Pro i^)'^ + Pmo ijf + PA late-time 

where ao is the current scale factor, pro, Pmo are p\ are 
the current energy densities corresponding to radiation, 
matter and the cosmological constant, respectively. On 
the other hand, at the radiation dominated epoch the 
energy densities (j2.12p are equal implying 



B = (proao)^'^" . 



(2.13) 



Finally, we notice that A is related to the scale of infiation 
(see Eq. (|23)) ). 



III. THE INFLATIONARY DYNAMICS 

The inflationary dynamics of the model presented in 
section|TT]is better described through a minimally coupled 
scalar field, 4>, with a potential, V^(0), whose energy 
density and pressure read 



2a2 



y(0) 



I0? 



y(0). (3.1) 



In the previous equations the prime stands for derivative 
respect to the conformal time. 




0.2 0.4 0.6 0.8 1 1.2 1.4 

+ «|0 



as ptji = p and = P. The last two equalities imply ^ 

coshT+° ^k|1 -t- aj (/)^ 
+2 cosh^T^ (^K|H-a|0) , (3.2) 



where 



(3.3) 



In addition, can also be determined analytically in 
terms of the scale factor 



(a)=± 



1 



K 1 + a\ 



rarcsinh 




-2(l+a) 



(3.4) 



For simplicity, we will restrict to the solution p.4p with -|- 
as our results does not depend on which of the two signs 
we choose. The universe starts its infiationary phase in 
a quasi-dS state where the scalar field is sitting close 
the top of the potential; i.e. at ^ (see Fig. [T|). 
Then it starts rolling down the potential until it exits 
the inflationary era when ipifli,) = 0*; i.e. 



ln(l + ^/2) 
k|1 + al 



(3.5) 



Finally, the radiation dominated phase starts for large 
values of the scalar field, i.e. 0* ^ (f). For latter 
convenience, we derive as well in terms of the scalar 
field 



= Vo coshi+° (k|1 -I- Q!|0) . 



(3.6) 



The number of e-folds of expansion since a given mode 
k exits the horizon (fc — aH) during the infiationary era, 
ai (j) — (pci until the end of infiation can be approximated 
by 6 



(3.7) 



Substituting the scalar field potential (13. 2p on the 
previous expression, we obtain 



1 



2(l + a)(2a- 1) 
1 - 4 a + cosh (2k| 1 -|- a 



4(1 - a) 




FIG. 1: Plot of the potential (|3.2p as a function of the scalar 
(j) given in Eq. (|3.4p. 

In a Friedmann-Lemaitre-Robertson- Walker (FLRW) 
universe, the scalar field, (/>, is equivalent to the modified 
Chaplygin gas introduced in Eqs. (|2.1I) and (|2.2p . as long 



^ A similar potential for a different equation of state was suggested 
in Ref. 25] . 

" Here, we are making use of a slow roll approximation |26|I : i.e. 
values of the scalar field such that e, |t;| <IC 1 where e and r) are 
the standard slow-roll parameters [2al (see also Eq. 114. 4| |). 
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IV. OBSERVATIONAL CONSTRAINTS 

We constrain our model through 
CMBR/LSS measurements, using the vahies 
P3(fcc) = (2.45 ±0.23) X 10-9 and = 1.0 ± 0.1 

for the power spectrum and spectral index of the 
scalar perturbations, respectively, evaluated at 
the pivot wave number kc = 0.05 Mpc"^, that 
corresponds to CMBR/LSS scales that exited the 
horizon during inflation (26j . To first order in the 
slow-roll approximation, the amplitude of the power 
spectrum for density perturbations is given by [26l | 



Ps{k) 



127r2 Vl' 



(4.1) 



k=aH 



where = Notice that Pg is proportional to Vq 

and therefore the power spectrum is extremely useful to 
constrain the energy scale of infiation (cf. Eqs. (|2.3I) . 
(|3.2p and (|3.3I) '). We will follow two different methods to 
constrain Vq: 

• Method 1: As a first step to constrain our model, 
we use the bounds on Nc, i.e. the number of e- folds 
since a given mode exits the horizon at 0c until 
the end of inflation. We consider as a conservative 
range 47 < A^c < 62 [11]. Then for a given 
Nc, on the allowed range, and for a fixed a, we 
obtain the corresponding value (j)c (cf. Eq. p.Sp V 
Finally, using the pair (a, (j)c) in Eq. (|4.ip with 
Ps{kc) = (2.45 ± 0.23) X 10"^, we deduce the scale 
of inflation; i.e. Vq. This is precisely, what is shown 
on the dashed-dotted and dotted lines of Fig. [51 

• Method 2: A tighter constraint on Vo is obtained 
by evaluating Pg at the pivot scale kc when the 
mode exits the horizon; i.e. kc = aH, during the 
inflationary era. Using Eqs. (j2.6p . (|3.4p and (|3.6I) . 
this scale reads 

27r 1/4, ,1/4 

kc = -^'^Pro oo X 



X ( cosh(K|l + a\4'c) coth(K|l + a\ 



.(4.2) 



Finally, by combining the previous equation and 
Eq. (|4.ip . we deduce the corresponding scale Vq for 
a given parameter a. The results are shown on the 
solid line in Fig. O 



We next consider the tensor to scalar ratio 



Ps{k) 



= 16e, 



(4.3) 



where Pt{k) is the tensorial power spectrum [2^]. We 
analyse the variation of this ratio against the spectral 
index of the scalar perturbations = 1 — 6 e + 2 r]. In 
the previous expressions the slow-roll parameters read 



r?(a,0c) = A^. (4.4) 

V 
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-1.03 -1.02 
a 



-1.01 
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FIG. 2: Variation of the inflationary scale versus a. Using 
the first method, i.e. the number of e-folds and Ps, we obtain 
the dashed-dotted and dashed lines corresponding to A*'c = 47 
and A'^c ~ 62, respectively. Using the second method; i.e. Ps 
at the pivot scale kc — 0.05 Mpc"^, we obtain the sohd hue. 



We follow the two approaches enumerated before to 
evaluate r and Hs . Our results are shown on the — r 
parameter space in Fig. [3l We see that the spectrum 
is slightly more red than preferred by the observation. 
However, Fig. [3] shows that r < 1 is in agreement with 
the observation 26]. Our best fit value corresponds to 
a — —1.024; i.e. an energy scale for inflation about 
10^^ GeV. It is worthy to notice that for a standard 
Chaplygin gas accommodated to the late-time expansion 
of the universe, negative values of a are allowed by Type 
la Supernova'' [27| . 



0.01 



0.001 



0.0001 




FIG. 3: Constraints on the Us—r parameter space. Using the 
first method; i.e. i.e. the number of e-folds and Ps, we obtain 
the dashed-dotted and dashed lines corresponding to A^c = 47 
and A'^c ~ 62, respectively. Using the second method; i.e. Ps 
at the pivot scale kc — 0.05 Mpc"^, we obtain the solid hue. 
The black lines are for a constant a-value, parameterised by 
an increase in the number of e-folds. 



We thank O. Bertolami for pointing out this to us. 
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GRAVITATIONAL WAVE SPECTRUM 



expansion that takes the form 28 1 



The stochastic gravitational wave (GW) spectrum, 
generated during the expansion of the universe, can be 
calculated using the method of continuous Bogoliubov 
coefficients, using the methods first developped by L. 
Parker fl2l and A. A. Starobinsky [l^, and later 
applied in (l5l [28|. Indeed, as the universe evolves, the 
annihilation and creation operators change with time. 

We can write these operators in terms of time- 
fixed annihilation and creation operators through a 
Bogoliubov transformation, involving the Bogoliubov 
coefficients a and S, which must satisfy the constraint 
Iap-|/3P = 1 EEl- 

It can be shown that the dimensionless relative 
logarithmic energy spectrum of the GWs, ricw, at the 
present time 770, is related with /3 [28{ : 



GW 



1 dp, 



GW 



-uj^P\ (5.1) 



where pcw is the energy density of GWs and cj 
the respective angular frequency; pc and H are the 
critical density of the universe and Hubble parameter, 
respectively, evaluated at the present time. 

Therefore, the evolution of the Bogoliubov coeSicient 
/3, since early times, gives us the power spectrum at the 
present time. The same coefficient can be written in 
terms of two continuous functions of time, X and in 
the form — [X — y)^/4. These new functions are 
determined by the set of equations [ll] 



X' = -ikY, 
Y' = 



fc V a 



X. 



(5.2) 
(5.3) 



The Friedmann equation (j2.6|) and the conservation 
of the energy momentum tensor imply a" /a cx p — Zp. 
Therefore, the power spectrum depends on the matter 
content of the universe (see for example Eq. (15.31) '). 
Consequently, the GWs power spectrum can provide us 
with important information about the universe evolution 
and its matter content, as long as the latter leaves a 
significant imprint on the spectrum. 

As our model describes a transition between the 
inflationary period and the radiation-dominated era, 
the integration must begin at the inflationary period, 
where we have a constant energy density. Hence, it is 
necessary to specify, at this inflationary period, the initial 
conditions for X{r}) and Y{r]). 

To solve the set of differential equations (|5.2p and (|5.3p , 
for our model, we must use numerical methods, but, since 
it admits an exact analytical solution for the case of a dS 
universe [28j . i.e. similar to our model at early-time (see 
Eq. (j2.3p ). we can specify appropriate initial conditions 
for X{rj) and Y{7]), at an appropriate scale of inflation. 
In this case, Eqs. (|5.2I) and (j5.3|) yield a solution for a dS 



= ( 1 + ^^ I e 



fc2 



(5.4) 
(5.5) 



where at is a scale factor during the dS phase. We use 
these expressions as initial conditions for the numerical 
integration of Eqs. ([0|) - ((0)) . 

The other piece of information we need to solve 
Eqs. (I5.2p - (|5.3|) is the whole evolution of the universe: 



a 
a 



early-time 



(5.6) 



late-time 



where the transition between the two epochs is 
determined by the moment when o"/o (for both periods) 
is approximately equal at a = aint . We use the values for 
the energy densities parameters given by WMAP5 '^^i 
a„ = 0.291 ± 0.014 and f^A = 0.709 ± 0.014, through the 
usual definition $7 = -Sr^p, where i?o — 71.3 km/s/Mpc. 

For a given mode k, the integration is done in terms 
of the scale factor a. We start at a < a^, just before 
the mode exits the horizon; i.e. k < ku where k^ = 
aH, until a ^ af, when the mode reenters the horizon; 
i.e. k < kn (c.f. Fig. |4|). The maximum value of k^ 
corresponds to the maximum height of the potential a" /a 
as the gravitons are produced only below the barrier a" /a 
(see Fig. g]) 
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FIG. 4: A scheme of the integration method used in the 
numerical simulations for a given mode k. The solid line refers 



to 



and the dashed-dotted one to kjj = {aH) 



To ensure the accuracy of the numerical results are 
not a numerical artifact we use the condition on the 
Bogoliubov coefficients — |/3p = l, which in terms of 
the new variables X, Y reads 2XY = 1 , as a constraint 
equation of the numerical solution. 
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VI. NUMERICAL SIMULATIONS 

The GW spectra resulting from the numerical 
integrations for the two different methods presented in 
section IV to constrain the model observationally are 
shown in Figs [5] and [6l In the former method, we use 
as a constraint equation the number of e-folds, N^, and 
in the latter the mode function, kc. Therefore, as one 
would expect, the variation with the values of a is the 
same, independently of the method used to calculate the 
inflationary scale which gives us the constant A of our 
model (see Eq. ([331) and Fig. ^. As we can see m 
Figs [S] and ini the spectrum changes with a, in particular 
it suffers simultaneously a vertical displacement and a 
variation in the high frequency region. 




f (Hz) 

FIG. 5: Gravitational wave spectra for the values a = 
— 1.05 (black), —1.04 (gray), and — 1.03 (light gray), along 
with Nc = 47 (solid) and Nc = 62 (dashed). The 
spectrum suffers a downward shift with an increase in |q:|, 
and, simultaneously, an increase in the high frequency region. 

Furthermore, to understand the vertical displacement 
of the plateau in the intermediate region of the spectrum, 
it is necessary to consider the inflationary scale as a 
function of a, plotted in Fig. [S] The same figure shows 
that, for a constant, an increase in the number of e-folds 
from Nc = 47 (dashed-dotted line) to = 62 (dotted 
line) results in a decrease in the inflationary scale Vq. 
Therefore, the decrease of the plateau height* , which can 
be observed in Fig. [S] by the difference between the solid 
lines {Nc = 47) and the dashed lines {Nc = 62) for a 
constant a (a given type of line) , is easily interpreted as 
a decrease in the inflationary scale. 

On the other hand, having showed that the vertical 
displacement of the spectra depends only on the value 
of the inflationary scale, we immediately conclude that 
the increase in the high frequency region, when a moves 
away from —1, is a function of the parameter a only, as 
we can see for example in Fig. [B] 



An early dS phase results in a flat plateau in the spectra. The 
height of the plateau is an increasing function of the Hubble rate 
of the dS expansion; i.e. the energy scale of the dS inflationary 
era. 



The results presented above show that, in the 
gravitational-wave energy spectra for the unified 
inflationary and radiation eras consistent with 
CMBR/LSS constraints, and in particular the one 
by means of the mode function in Fig. |6l present a 
variation in the spectrum that is a combination of the 
last effects: a vertical shift as a result from the decrease 
in the inflationary scale; and an increase in the high 
frequency region that results from a increase in the 
|a|-value. 

The last variation is the more important one because 
it can give us a new insight about the relation between 
the frequency limits of the GW spectrum and the model 
behind the inflationary and the radiation eras. 




10-18 10-14 10-10 10-6 0.01 100 



f (Hz) 

FIG. 6: Gravitational wave spectra for the values a = 
-1.05 (black), -1.04 (gray), and - 1.03 (light gray) . The 
spectrum suffers a downward shift with an increase in |q:|, 
and, simultaneously, an increase in the high frequency region. 



Although the spectra presented on this section are 
smooth, reflecting the smoothness of the transition 
between the inflationary and the radiation dominated 
periods, characteristic of the present model, this is not 
always the case. For instances, when the system goes 
through a complex transition period, with changes in the 
equation of state, the structure of the spectrum is equally 
complex, particularly in the hundreds of MHz and GHz 
regime, allowing us to read a lot of information about the 
properties of the transition (see ref. [28]). 

What about the reheating of the universe in our model? 
such a process can be modelled by coupling the scalar 
field to a radiation fiuid, where the decay rate of the 
scalar field into the radiation fluid is governed by a 
phenomenological parameter. In addition, the tail of 
the potential p.2p is substituted by a potential with a 
vanishing minimum. For the modified model the trend 
and the level of the spectra of the GW would not be 
much different to the one presented on this section. 
Indeed, it is still the parameter a that characterises the 
spectra. In addition, the inclusion of a reheating process 
is computationally very much time consuming. 
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VII. CONCLUSIONS 

In this work, we investigate the production of GWs 
in a new modified generalised Chaplygin gas for the 
early universe (cf. Sees [TT] and IIIip . Through recent 
measurements of CMBR/LSS, the parameters of the 
model have been constrained. 

We have used the method of the continuous Bogoliubov 
coefficients to calculate the GW energy spectrum for 
different scales of energy and values of the a-parameter 
of the model. 

Besides the fact that the model suffers a small 
deviation from the preferred values of the spectral 
index n^, the obtained spectra reveal a consistent 
picture corresponding to a smooth transition between 
the inflationary and the radiation dominated epochs of 
the universe, and constitute a significant imprint of this 
modified generalised Chaplygin gas. 

Finally, the strong variation of the high frequency 
range, is directly related with the decrease in the 
maximum of the potential a!' ja (see Eq. (|5.3p ) for |a| 
approaching 1. This feature implies strong limits to the 
maximum frequency allowed in our model and which 
will be within the reach of future gravitational-waves 



detectors like BBO and DECIGO for the KHz range 
of frequencies. In fact, for the most consistent values of 
the a-parameter, the spectrum shows a frequency as low 
as in the Hz region. 

Concluding, these results show that, for this model, 
the Hz-KHz frequency range comes directly from the 
transition between the inflationary regime and the 
radiation era, and addresses important issues about the 
limits of the GW energy spectrum. 

Last but not least, we would like to highlight that 
one of the merits of the model we have presented for 
the transition from the inflationary era to the radiation 
dominated epoch is its relative simplicity. 
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